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1.  Introduction 

The  Back-EMF  theory  of  Onset  has  been  refined  and  expanded  to  enable  prediction  of  erosion 
rates  for  steady  state,  self-field  MPD  thrusters[l],  This  theory  is  capable  of  explaining  both 
the  observed  oscillations  as  well  as  the  increased  erosion  at  Onset.  This  represents  the  first 
time  that  the  plasma  discharge  and  electrode  processes  have  been  coupled.  Recent 
experiments  at  the  Institut  fur  Rahmfahrtsysteme  at  Stuttgart  reveal  cathode  damage  at 
precisely  the  same  location  as  that  predicted  by  this  theory[2].  The  research  described  in  this 
report  provides  a  summary  of  the  accomplishments  and  progress  made  during  the  three  years 
and  specifically  during  the  final  year  under  this  grant  AFOSR-87-0360. 

2.  Research  Objectives 


The  purpose  of  this  research  is  to  understand  and  quantify  the  mechanisms  responsible  for 
evaporative  erosion  in  steady  state  magnetoplasmadynamic  (MPD)  thrusters.  This  is  an 
important  step  in  being  able  to  predict  thruster  lifetimes.  Erosion  processes  depend  on  a 
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complex  coupling  between  plasma  discharge  characteristics,  plasma-wall  interactions,  and 
electrode  phenomena.  In  particular,  erosion  rates  depend  on  whether  the  current  conduction  is 
through  localized  spots,  or  via  a  diffuse  (distributed)  mode.  Spots  are  detrimenta,'  to  the 
electrode  material  because  of  their  high  erosion  rates.  Therefore,  it  is  important  to  understand 
how  and  under  what  conditions  they  may  be  formed,  and  exactly  when  diffuse  mode  behavior 
ends.  Much  of  the  focus  of  this  research  is  on  the  cathode.  The  following  overall  objectives 
were  set  for  this  research  effort:  (1)  to  identify  the  limits  of  steady  state  diffuse  mode 
electrode  operation,  (2)  to  identify  the  mechanisms  responsible  for  limiting  the  diffuse  mode, 
(3)  to  model  the  thermal  response  of  the  cathode  at  steady  state,  thereby  enabling  the 
prediction  of  surface  temperatures  and  erosion  rates  due  to  evaporation,  and  (4)  to  enable 
prediction  of  the  temporal  thermal  behavior  of  the  cathode.  Although  no  transient  model  of 
the  bulk  plasma  flow  in  the  MPD  thruster  exists  presently,  it  should  be  forthcoming  in  the 
near  future. 


3.  Status  of  the  Research  Effort 

This  section  describes  the  highlights  of  the  research  during  the  final  year  of  this  grant.  A 
more  detailed  discussion  can  be  found  in  the  appendices  as  well  as  in  the  references  at  the 
end  of  this  report. 

The  Back-EMF  theory  was  first  formulated  within  the  context  of  frozen,  fully  ionized  channel 
flow  in  the  MPD  thruster[3].  It  was  found  that  MPD  channel  flow  is  parametrized  by  a  non- 
dimensional  parameter  S*  =  B*1/\i0Fa* ,  called  the  Magnetic  Force  Number.  S*  is  defined 
with  respect  to  quantities  evaluated  at  the  magnetogasdynamic  sonic  point  (i.e.  where  M*=l), 
where  B*  is  the  magnetic  induction,  (i0  is  the  permeability  of  free  space,  F  is  the  mass  flux, 
and  a*  is  the  frozen  speed  of  sound  at  the  sonic  point.  S*  may  be  interpreted  as  the  ratio  of 
magnetic  pressure  to  kinetic  energy  density.  Careful  analysis  of  MPD  channel  flow  revealed 
that: 

(1)  S*  is  related  to  the  experimentally  measured  Onset  parameter  J2/m, 

(2)  fundamental  limits  exist  on  S*  and  hence  on  J2/m  for  supersonic  flow  to  exit  the  MPD 
channel,  and 


(3)  Onset  limits  could  be  derived  from  first  principles,  in  agreement  with  the  experiments  of 
Malliaris  et.  al.[4j. 
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Incorporation  of  finite  rate  ionization  and  recombination  did  not  alter  these  results,  but  an 
explicit  analytical  closed  form  expression  for  the  limits  on  S*  could  not  be  found[51. 
However,  the  numerically  calculated  Onset  limit  predicted  by  this  theory  agreed  with  quasi¬ 
steady  experiments  in  straight  coaxial  thrusters[6],  The  Onset  limit  predicted  by  this  theory 
was  related  to  excessive  values  of  the  Back-EMF  which  occur  at  high  values  of  the  total 
current  for  a  given  mass  flow  and  geometry.  This  steady,  quasi-one  dimensional  flow  model 
provides  variations  of  the  important  quantities  (number  densities,  current  density,  temperature, 
velocity,  ionization  fraction,  and  magnetic  induction)  versus  distance  along  the  channel.  This 
is  dene  for  a  given  mass  flow  rate,  geometry,  and  total  current. 

The  plasma  flow  and  electrode  processes  in  the  MPD  thruster  are  intimately  coupled.  A  first 
attempt  at  treating  this  complex  coupling  consisted  of  using  the  axial  profiles  generated  by  the 
quasi-one  dimensional  model  as  boundary  conditions  for  a  two- temperature  model  of  the 
electrode-adjacent  boundary  layer  flow[7].  The  boundary  layer  model  provided  distributions 
of  heavy  particle  temperature,  ionization  fraction,  and  number  densities  in  the  direction 
transverse  (i.e.  radial)  to  the  electrodes.  Two  important  outcomes  of  this  boundary  layer 
theory  were  (a)  that  the  heavy  particle  conductive  heat  flux  was  negligible  compared  to 
charged  particle  bombardment  of  the  electrodes,  and  (b)  that  the  highly  viscous  MPD 
boundary  layer  flow  produces  a  maximum  in  the  transverse  heavy  particle  temperature 
distribution  somewhere  within  the  boundary  layer  well  away  from  the  core  flow  and  the 
electrode  surfaces.  The  boundary  layer  theory  provided  important  boundary  conditions  for 
analyses  of  the  sheath  regions  near  the  cathode  and  anode[8,9].  Analyses  of  the  sheaths 
provided  the  final  set  of  boundary  conditions  for  thermal  analyses  of  the  electrodes[10,l  lj. 
These  simple  sheath  models  revealed  several  remarkable  phenomena: 

(1)  Stable,  diffuse  mode  behavior  would  transition  into  the  spot  mode  at  a  critical  value  of  the 
sheath  voltage  drop. 

(2)  This  diffuse  to  spot  transition  was  caused  by  a  thermal  runaway  due  to  excessive  electron 
bombardment. 

(3)  The  thermal  runaway  was  found  to  occur  in  the  low  current  density  regions  of  the  cathode 
and  the  high  current  density  regions  of  the  anode. 

(4)  As  the  total  current  to  the  thruster  is  increased  for  a  given  mass  flow  and  geometry,  the 
current  density  profile  along  the  thruster  length  reaches  a  minimum  (due  to  excessive  Back- 
EMF)  in  the  middle  while  reaching  local  maxima  at  the  inlet  and  exit. 
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(5)  The  regions  of  the  anode  closest  to  the  inlet  and  exit  of  the  thruster  are  thus  at  greatest 
risk  in  terms  of  erosion,  while  the  middle  region  of  the  cathode  is  susceptible  to  erosion  at  a 
critical  value  of  the  current. 

(6)  The  present  theory  also  predicts  that  unstable  voltage  and  electric  field  oscillations  may 
occur  at  the  high  current  density  regions  (i.e.  near  inlet  and  exit  regions)  of  the  cathode. 
These  oscillations  may  be  triggered  by  the  energetic  electrons  emitted  by  the  cathode  surface 
interacting  with  the  relatively  slower  plasma  electrons  (classical  two-stream  instability). 
Hence,  the  regions  of  greatest  erosion  are  predicted  not  necessarily  to  occur  in  the  regions  of 
greatest  oscillations. 


Recently,  there  were  two  reports  of  tests  conducted  at  Stuttgart  on  high  power,  steady  state 
MPD  thrusters[2].  These  tests  showed  that  in  the  two  different  geometries  tested(one  with  a 
flared  anode  and  the  other  with  a  straight  anode),  severe  damage  was  encountered  at  the 
respective  cathodes.  In  fact,  the  cathodes  appear  to  have  exploded  from  within  (see  Fig.  1). 
This  damage  occured  at  6500  A  for  the  flared  anode  (DTI  thruster)  and  at  8000  A  for  the 
straight  coaxial  ZT1  thruster.  The  Back-EMF  theory  of  Onset  not  only  correctly  predicted 
these  current  limits  but  also  predicted  the  temperature  profiles  shown  in  Fig.  2.  As  expected 
the  centerline  cathode  temperatures  far  exceed  the  cathode  surface  temperatures  at  the  Onset 
limit.  Furthermore  at  the  predicted  Onset  limit,  the  interior  temperature  exceeds  the  melting 
temperature.  When  this  occurs,  pressure  builds  up  inside  the  cathode  driving  crack 
propagation  and  promoting  fissures.  Although  the  Back-EMF  theory  correctly  Dredicts  this 
Onset  and  damage  limit,  it  cannot  describe  the  events  after  the  damage  and  melting  have 
begun  to  occur. 

In  summary,  the  simple  models  developed  thus  far  enable  the  designer  to  use  quick  tools  for 
performance  studies  and  design  evaluations.  Efforts  are  underway  to  develop  quantitatively 
accurate  predictions  of  the  flow  field  in  straight  coaxial  thrusters  (2D  axisymmetric 
simulations).  Accurate  models  are  still  needed  of  the  electrode  pre-sheath  regions  where  the 
prediction  of  ionization  fraction  is  important.  Specifically,  research  is  underway  to  examine 
the  combined  effects  of  trapping  of  resonance  radiation  and  ionization  of  excited  species  by 
electrons  emitted  from  the  electrodes. 
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Fig.l:Scht.  matic  of  the  type  of  damage  sustained  by  the  DT2  and  ZT1  thrusters  after  sustained 
steady  state  operation  at  6500  A  and  8000  A  respectively.  The  total  argon  propellant 
flow  reported  was  2  g/s.  For  further  details,  please  see  reference  [il]. 


Fig.2:Shown  here  are  calculated  surface  and  centerline  temperature  profiles  versus  distance 
along  the  cathode  for  the  Stuttgart  ZT1  thruster.  The  horizontal  line  at  3680  K  denotes 
the  melting  temperature  of  tungsten.  Temperatures  above  this  value  cannot  be  taken 
seriously  from  the  theoretical  predictions,  since  the  theory  does  not  account  for 
phenomena  once  melting  has  occurred.  However,  as  the  model  suggests  and  as  shown  in 
this  figure,  the  centerline  temperature  has  exceeded  the  surface  temperature  in  the  select 
region  thereby  explaining  the  exloded  cathode  observed  in  the  Stuttgart  t^st. 
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Abstract 

Electrode  erosion  in  MPD  thrusters  is  significantly  dependent  on  whether  the  mode  of 
operation  is  diffuse  (i.e.  the  net  current  density  is  distributed  over  the  electrode 
surface)  or  via  spots.  Since  spot  mode  erosion  rates  are  typically  much  higher  than  the 
difuse  mode  erosion  rates,  the  latter  is  more  desirable  from  the  standpoint  of  sustained 
steady  state  operation.  In  this  paper,  the  steady  state  thermal  response  of  the  cathode 
in  an  MPD  thruster  is  analyzed  in  order  to  determine  the  limits  of  diffuse  mode 
operation.  An  energy  balance  at  the  cathode  surface  along  with  current  conservation, 
serve  to  simultaneously  determine  the  surface  temperature  and  the  total  sheath  voltage 
drop  at  a  given  location  on  the  electrode.  It  is  found  that  two  different  limits  to  steady 
state  diffuse  mode  operation  exist.  One  limit  corresponds  to  an  unsteady  thermal 
runaway  caused  by  excessive  electron  bombardment,  in  which  regenerative  heating 
leads  to  local  melting  of  the  electrode  surface.  The  other  limit  corresponds  to  an 
operating  regime  where  a  steady  state  cannot  exist  because  of  increased  ion 
bombardment.  In  this  regime,  a  small  increase  in  the  net  plasma  current  density 
results  n  a  large  sheath  voltage  drop.  This  causes  an  increase  in  local  surface  heating 
through  ion  bombardment,  until  a  new  state  is  reached  which  is  inherently  unsteady 
and  accompanied  by  plasma  oscillations.  These  two  limits  found  here  may  correspond 
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to  transitions  from  diffuse  to  pre-spot  modes  of  electrode  operation.  These  operating 
limits  are  also  found  to  be  strongly  influenced  by  electrode  and  discharge  parameters, 
as  well  as  by  external  cooling. 


Nomenclature 

A  Thermionic  emission  constant 

E  Electric  field 

e  Electronic  charge 

h  Heat  transfer  coefficient  roughly 

representing  overall  heat  transfer  to  an 
external  coolant  which  is  at  a 
temperature  Tcool 

k  Boltzmann’s  constant 

L  Cathode  length 

Q  Net  heat  per  unit  area  per  unit  time  into 

the  cathode  surface,  at  a  given  axial 
location 

T  Local  cathode  surface  temperature 

V  Voltage  or  potential 

y  Coordinate  in  the  direction  normal  to  the 

cathode  surface 

£,•  Ionization  potential  of  the  propellant  gas 

<}>  Electrode  work  function 

Ce  Absolute  velocity  of  an  electron 

Component  of  the  absolute  velocity  of 
an  electron  in  the  y-direction,  i.e.  in  the 
direction  normal  to  the  surface 

Ec  Magnitude  of  the  y-component  of  the 

electric  field  at  the  cathode  surface 

Ew  y-component  of  the  electric  field  at  the 

sheath-edge 

fe  Electron  velocity  distribution  function 

jg  Current  density  of  plasma  electrons 

jg  Surface  electron  emission  current  density 
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ji  Cuirent  density  of  plasma  ions 

y'oo  Net  plasma  current  density  at  a  given 
axial  location 

Kt  Thermal  conductivity  of  cathode 

material 

me  Mass  of  an  electron 

mt  Mass  of  an  ion 

ne  Electron  number  density 

nt  Ion  number  density 

n„  Charged  particle  number  density  at  the 
edge  of  the  collisionless  sheath,  at  a 

given  axial  location 

> 

T cool  Temperature  of  an  external  coolant 

T „  Plasma  electron  temperature  at  a  given 
axial  location 

Vc  Value  of  the  voltage  at  the  cathode 

surface  at  a  given  axial  location,  or  the 
total  voltage  drop  across  the  sheath 
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1.  Introduction 

Electrode  erosion  is  of  primary  importance  for  the  prediction  of  lifetimes  of 
magnetoplasmadynamic  (MPD)  thrusters.  Erosion  processes  depend  on  a  complex 
coupling  between  plasma  discharge  characteristics,  plasma-wall  interactions,  and 
electrode  phenomena.  In  particular,  erosion  rates  depend  on  whether  the  current 
conduction  is  through  spots,  or  via  a  diffuse  (distributed)  mode.  The  diffuse  mode  is 
characterized  by  distributed  current  emission  over  the  electrode  surface,  surface 
temperatures  well  below  the  material’s  melting  temperature,  and  non-negligible  plasma 
ion  current.  By  contrast,  spots  are  characterized  by  lo^lly  constricted  or  filamentary 
current  conduction,  surface  temperatures  at  or  above  the  material’s  melting 
temperature,  and  strong  temperature-field  emission.  Spots  are  detrimental  to  the 
electrode  material  because  of  their  high  erosion  rates[l].  Therefore,  it  is  important  to 
understand  how  and  under  what  conditions  they  may  be  formed,  and  exactly  when 
diffuse  mode  behaviour  ends.  Current  understanding  of  this  transition  from  diffuse  to 
spot  behaviour  is  at  best,  rather  poor.  In  this  paper,  the  steady  state  diffuse  mode  at 
the  cathode  is  examined  theoretically.  Limits  on  thruster  operation  in  this  mode  are 
found,  and  an  explanation  for  a  possible  transition  to  spot  formation  is  given. 
Additionally,  the  analysis  yields  a  means  of  estimating  cathode  surface  temperature 
which  can  be  used  to  predict  evaporative  erosion  rates  for  the  cathode. 

Earlier  research  on  cathode  processes  has  focused  mainly  on  spot  behaviour  and  spot 
characteristics[2-10].  Although  these  works  reveal  the  intricate  and  complex 
phenomena  occuring  in  a  cathode  spot,  little  or  no  information  is  provided  as  to  how 
spots  are  formed  in  the  first  place.  Some  authors  have  attempted  to  describe  transition 
from  diffuse  to  spot  modes.  Moizhes  and  Rybakov[ll]  have  found  a  negative-slope 
region  in  the  emission  current  voltage  characteristic,  which  they  attribute  to  transition 
to  a  spot  mode.  These  authors  attribute  the  transition  to  a  thermal  instability,  whose 
physical  meaning  or  origin  is  totally  unclear.  A  thermal  runaway  mechanism  in  spots 
due  to  a  positive  feedback  between  Joule  heating  and  a  temperature  dependent 
electrical  conductivity,  has  been  proposed  by  Hantzsche[12].  In  the  same  paper 
however,  the  author  shows  that  thermionic  emission  cools  the  surface  and  prevents  this 
thermal  runaway.  It  is  important  to  note  that  Hantzsche  assumes  the  sheath  voltage 
drop  to  be  constant  and  given. 

By  contrast  to  the  cathode,  transition  at  the  anode  appears  to  have  been  studied  more 
extensively[13-19].  These  works  attribute  spot  formation  to  an  initial  local  surface 
meltdown  arising  primarily  from  Joule  heating  at  high  current  densities.  However, 
another  thermal  runaway  mechanism  involving  the  much  ignored  sheath  exists[20].  In 
this  theory,  electron  emission  from  the  electrode  surface  decreases  the  sheath  potential 
drop  resulting  in  increased  electron  bombardment  from  the  plasma.  This  results  in  a 
positive  feedback  between  increased  emission  and  heating  due  to  electron 
bombardment,  causing  a  thermal  runaway.  This  may  occur  when  the  anode  sheath 
voltage  difference  (defined  here  as  anode  potential  minus  the  adjacent  plasma 
potential)  is  negative  and  increasing.  This  thermal  runaway  mechanism  arises  before 
the  well  known  anode  sheath  reversal[21],  and  may  explain  formation  of  anode  spots. 
In  this  paper,  the  same  mechanism  is  shown  to  be  also  operative  at  the  cathode. 
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A  simple  model  of  the  cathode  surface  and  its  adjacent  sheath  is  discussed  in  section 
2.  Results  for  a  purely  thermionically  emitting  cathode  under  a  range  of  conditions 
corresponding  to  MPD  thruster  operation,  are  given  in  section  3.  The  physical 
meanings  of  the  limits  on  steady  state  diffuse  mode  operation  are  discussed  in  section 
4,  and  the  effects  of  field-enhanced  and  temperature-field  emission  are  discussed  in 
section  5.  Finally,  a  summary  of  this  work  along  with  its  conclusions  are  presented  in 
section  6. 


2.  A  Simple  Model 


This  section  will  focus  on  a  thermal  model  of  the  cathode  surface  at  steady  state.  The 
cathode  surface  is  subjected  to  both  heating  and  cooling  (Fig.  1).  Charged  particles 
(i.e.  electrons  and  ions)  from  the  plasma  bombard  the  surface  and  consequently  heat  it 
The  ions  (which  we  will  assume  to  be  singly  charged)  also  recombine  with  electrons  at 
the  electrode,  thereby  heating  the  surface  even  further.  Additional  heating  due  to 
radiation  from  the  plasma  is  also  present  and  expected  to  be  important[22],  but  will  be 
neglected  here.  The  electrode  may  also  emit  electrons  which  can  result  in  heating 
(during  field  emission)  or  in  cooling  (during  thermionic  emission).  During  steady  state 
operation,  it  is  likely  that  the  electrode  is  cooled  extemally[l,23].  This  is  usually  done 
by  cooling  the  cathode  at  its  base.  Finally,  the  cathode  surface  may  radiate  away  its 
heat  With  these  considerations,  we  may  proceed  to  write  the  following  for  the  net 
heat  per  unit  area  and  per  unit  time  entering  the  cathode  surface: 


Q  =Je 


2  kT„ 

(j>  + - 

e 


+  ji(Vc+Ei-$)  -  J'e 


-h(T-Tcool) 


(1) 


where  Q  is  the  net  heat  entering  the  surface  per  unit  area  and  per  unit  time,  je  is  the 
electron  current  density  from  the  plasma  into  the  cathode,  jt  is  the  ion  current  density 
from  the  plasma  into  the  cathode,  jE  is  the  current  density  of  emitted  elections  from 
the  surface,  e  is  the  charge  of  an  electron,  <{>  is  the  electrode  work  function,  T „  is  the 
temperature  of  the  plasma  electrons  at  the  edge  of  the  sheath,  T  is  the  surface 
temperature,  Vc  is  the  total  voltage  drop  across  the  sheath,  £,•  is  the  ionization  potential 
of  the  gas  (propellant),  h  is  an  overall  heat  transfer  coefficient,  and  is  the 
temperature  of  an  external  coolant.  The  quantity  hiJ-T^i)  represents  a  very  rough 
model  of  the  overall  cooling,  h  can  be  considered  as  equivalent  to  the  electrode 
material  thermal  conductivity  per  unit  length  (i.e.  KtIL).  At  steady  state,  Q  must  be  0 
and  equation  (1)  determines  the  surface  temperature  T.  Radiative  heat  transfer  to  the 
anode  can  be  shown  to  be  negligible,  and  is  therefore  ignored  here. 


Consider  next,  the  sheath  region  immediately  adjacent  to  the  electrode  surface.  This 
region  is  typically  of  the  order  of  tens  of  Debye  lengths.  The  electron  mean  free  path 
however,  is  orders  of  magnitude  larger[22].  Consequently,  the  sheath  may  be  treated  as 
a  collisionless  sheath.  The  electrons  entering  the  sheath  from  the  plasma  are  then 
described  by  the  collisionless  steady  Boltzmann  equation,  which  in  one-dimension  is: 


r  +  eE 


=  0 


(2) 
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where  fe  is  the  electron  velocity  distribution  function,  is  the  component  of  the 
absolute  velocity  of  the  electrons  in  the  direction  transverse  to  the  electrode,  y  is  the 
transverse  coordinate,  me  is  the  mass  of  an  electron,  and  E  is  the  electric  field.  Using 
the  following  transformations. 


%  = 


_  _eV_ 

2kT„  +  kT„ 

along  with  £  =  -  dV  Idy ,  yields 

dfe 


and 


t\  =  y 


=  0  ,  or  that  fe=fe(Z>)  only 


(3) 


(4) 


With  the  condition  that  the  electron  distribution  function  be  Maxwellian  at  the  sheath- 
edge  (where  V  is  taken  to  be  0),  we  obtain  the  following  exact  solution  for  the 
electron  distribution  function  in  the  sheath: 


fe  =  ne(mel2nkTJm  exp[-  meCe2/2kT„\  (5) 

where 

ne  =  exp[-  eVlkT„]  (6) 

with  being  the  electron  number  density  at  the  edge  of  the  sheath  where  V  is  taken 
to  be  0.  The  average  flux  of  electrons  may  then  be  easily  calculated  from  equations 
(5)  and  (6)  by  integrating  fe  over  dCa  and  dCa  from  -  «  to  «>,  and  over  dC^  from 
0  to  «.  This  flux  when  evaluated  at  the  wall  and  multiplied  by  e,  gives  the  electron 
current  density  at  the  wall: 

je  =  en„(kT  J2me  7t)1/2  exp[-eVclkTJ  (7) 

where  -Vc  is  the  potential  of  the  cathode  with  respect  to  the  sheath  edge  (where 
V  =  0),  and  je  is  taken  to  be  positive  in  the  direction  away  from  the  cathode. 


The  ions  entering  the  sheath  arc  assumed  to  be  monoenergetic.  The  Bohm  criterion 
for  a  stable  monotonic  sheath  is  therefore  satisfied [24].  This  condition  is  only  slightly 
altered  in  the  presence  of  electron  emission  from  the  surface[25].  More  exact  theories 
that  provide  the  velocity  distribution  function  of  the  ions  entering  the  sheath  from  the 
collisional  pre-sheath  exist[26-28].  However,  for  present  purposes  the  frequently 
employed  monoenergetic  assumption  simplifies  the  problem  a  great  deal.  We  may 
then  write  the  ion  current  density  at  the  surface  as: 

ji  =  en„(kr  Jmt)m  (8) 

where  7i  is  taken  to  be  positive  toward  the  cathode,  mi  is  the  mass  of  the  ions,  and 
quasi-neutrality  (n,  ~  ne  =  n„)  is  assumed  at  the  sheath  edge.  The  sheath  edge  is  a 
rather  arbitrary  definition  since  no  such  rigid  demarcation  exists  in  reality.  When  the 
ions  are  monoenergetic,  the  sheath  edge  has  a  clear  meaning.  However,  even  when 
the  ion  distribution  function  is  not  monoenergetic,  the  sheath  edge  can  be  defined  as 
the  location  where  the  velocity  of  the  ions  reaches  the  value  of  the  Bohm  velocity 
(kTJmi)1'2.  In  this  paper,  we  will  adopt  this  definition  and  take  the  datum  for  the 
sheath  potential  at  this  location  (i.e.  V  =  0  at  the  sheath  edge).  Now,  If  the  surface 
emission  current  density  denoted  by  y£  is  taken  to  be  positive  toward  the  surface  and 
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the  net  plasma  current  density  denoted  by  j ^  is  taken  to  be  positive  toward  the 
cathode,  we  may  write  overall  current  conservation  at  steady  state  as: 

j~  =  Je  +  ii  ~  Je  <9) 

where  in  general,  jE  is  a  function  of  the  local  surface  temperature  T  and  the  local 
electric  field  at  the  cathode  surface  Ec .  Determination  of  Ec  is  discussed  in  section  5. 
Equation  (9)  determines  the  total  cathode  sheath  voltage  drop  Vc.  Given  the 
parameters  n„,  T„,  4>,  h,  Tcooh  and  the  surface  electric  field  Ec,  equations  (1)  and 

(9)  represent  two  simultaneous  equations  in  the  two  unknowns,  T  and  Vc.  It  must  be 
pointed  out  that  Crawford  and  Cannara  mention  overall  current  conservation  in  their 
work[29].  However  their  regime  of  interest  was  at  very  low  densities  (10  m  ),  and 
hence  ion  and  electron  currents  from  the  plasma  were  negligible  in  comparison  with 
the  emission  current.  This  is  an  important  difference  between  their  work  and  the 
present  work. 


A  simple  solution  may  be  found  for  the  case  of  pure  thermionic  emission  (which  is 
likely  for  refractory  materials),  where  the  surface  emission  current  density  depends 
only  on  T.  For  this  case,  the  total  sheath  voltage  drop  may  be  explicitly  obtained 
from  equation  (9): 


h  +  Je  ~  Jo- 
Jr 


(10) 


where  jr  =  en„(kT J2nme)m.  Substituting  (10)  in  (1)  yields  the  following: 


Q 


2k  jE 
e 


2kT 

(T„-T)  +  y'(Vc  +e;+—=-) 


2  kT 

-  j~  (9  + - “)  -  h  ( T  -Tcool )  (11) 

e 

Since  at  steady  state  Q  must  be  zero,  equation  (11)  represents  an  implicit  equation  for 
the  local  surface  temperature  T,  where  Vc  is  given  by  (10). 


An  interesting  phenomenon  is  revealed  by  equation  (11).  The  emission  current  density 
which  under  thermionic  conditions  results  in  energy  transport  away  from  the  cathode, 
appears  as  an  input  heating  source.  This  can  be  understood  upon  close  examination  of 
equation  (9)  or  (10).  When  the  emission  current  density  increases  for  a  fixed  plasma 
number  density  and  plasma  current  density,  the  sheath  voltage  drop  must  decrease  in 
order  to  conserve  total  current.  This  results  in  an  increase  in  electron  bombardment, 
which  contributes  to  heating  the  surface.  Thus,  although  thermionic  emission  cools  the 
surface,  its  overall  effect  is  to  heat  the  surface  via  increased  electron  bombardment.  It 
is  therefore  important  to  include  the  effects  of  the  sheath  when  attempting  to  determine 
the  surface  temperature. 

A  simple  model  has  been  developed  in  this  section,  which  considers  a  thermal  balance 
of  the  cathode  surface  along  with  overall  current  conservation.  The  local  surface 
temperature  and  overall  sheath  voltage  drop  can  be  determined  from  this.  Results  for 
a  purely  thermionically  emitting  cathode  under  typical  MPD  conditions  are  presented 


next. 
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3.  Results  for  Pure  Thermionic  Emission 

The  governing  equations  presented  in  the  preceding  section  are  solved  in  this  section 
for  the  case  of  pure  thermionic  emission  (i.e.  no  field  enhanced  or  temperature-field 
emission).  Temperatures  and  sheath  voltage  drops  are  calculated  for  typical  MPD 
conditions. 

For  pure  thermionic  emission,  the  emission  current  density  is  given  by: 

jE  =  AT2  exp  [-e  §/kT  ]  (12) 

where  A  is  the  emission  coefficient  and  0  is  the  material  work  function.  The  sheath 
voltage  drop  Vc  can  then  be  written  in  terms  of  T  using  equation  (10).  Given  the 
parameters  n„,  jM  0,  h,  and  Tcoot,  equations  (10),  (11)  and  (12)  may  be 
combined  and  solved  for  T.  The  sheath  voltage  drop  can  then  be  determined  from 
(10).  The  emission  coefficient  may  vary  considerably  for  materials  with  oxide  layers. 
Both  A  and  0  may  vary  locally  on  a  given  surface,  so  that  it  is  important  to  consider 
the  sensitivity  of  final  results  to  variations  in  these  parameters.  This  is  addressed  later 
in  this  section. 

A  typical  variation  of  the  net  heat  into  the  cathode  (Q )  versus  the  surface  temperature, 
is  displayed  in  Fig.  2.  Two  intersections  with  the  horizontal  axis  are  found, 
representing  two  possible  steady  state  solutions.  The  first  (lower  temperature)  is  a 
stable  attractor,  while  the  second  (higher  temperature)  is  an  unstable  repeller.  This  can 
be  seen  quite  easily  by  perturbing  the  solutions  to  either  side,  and  determining  if  the 
initial  state  is  restored.  The  physical  meaning  of  the  stable  point  is  clear.  The 
incoming  energy  on  the  surface  is  exactly  balanced  by  the  outgoing  energy. 
Furthermore,  the  surface  can  maintain  this  temperature  for  an  indefinite  period  of  time 
even  for  small  fluctuations  in  the  discharge.  The  unstable  point  requires  some 
explanation.  Although  an  exact  balance  between  incoming  and  outgoing  energy  is 
possible  at  this  point,  the  surface  cannot  remain  in  this  state  indefinitely  and  is 
susceptible  to  change  due  to  infinitesimal  fluctuations.  Physically,  this  state  represents 
an  operating  regime  where  the  temperature  and  hence  the  emission  current  density  is 
high  enough  to  lower  the  sheath  voltage  drop  below  a  critical  value.  This  causes 
excessive  heating  due  to  increased  electron  bombardment,  resulting  in  a  temperature 
rise.  The  subsequent  increase  in  temperature  causes  an  increase  in  current  emission 
which  lowers  the  sheath  voltage  drop  further.  This  positive  feedback  process  then 
repeats  itself  until  the  surface  is  regeneratively  heated  upto  the  melting  temperature. 
This  effect  is  quantitatively  described  in  section  4. 

Results  from  the  present  calculation  are  shown  in  Figs.  3-6.  In  these  figures,  the 
steady  state  surface  temperatures  calculated  from  equation  (11)  are  shown  versus  the 
net  plasma  current  density  j„,  the  charged  particle  number  density  at  the  sheath  edge 
nM  the  electrode  work  function  0,  and  the  temperature  of  the  external  coolant  Tcool  • 
The  heat  transfer  coefficient  h  is  fixed  at  20  kW/m2/K,  the  plasma  electron 
temperature  is  taken  to  be  12000  K,  and  A  -  3  x  104  Aim2 IK1.  The  value  of  h  is 
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obtained  from  an  estimate  of  the  thermal  conductance  through  a  tungsten  cathode.  All 
four  plots  display  both  the  stable  (lower  temperature)  and  unstable  (higher  temperature) 
regimes.  Figures  3  and  4  display  the  variation  of  T  with  j„  and  respectively  for 
<))  =  2.63  V  and  Tcool  =  500  K.  From  Fig.  3,  it  can  be  seen  that  for  a  given  number 
density,  the  stable  diffuse  mode  is  limited  both  at  low  as  well  as  high  values  of  the 
current  density.  Similarly,  for  a  given  current  density  there  are  lower  and  upper  limits 
on  the  range  of  allowable  number  densities.  At  low  current  densities  («„  fixed)  and  at 
high  number  densities  (j „  fixed),  the  stable  and  unstable  solutions  are  seen  to  merge 
and  disappear.  Also,  at  the  higher  current  densities  (n„  fixed)  and  the  lower  number 
densities  (jm  fixed),  the  stable  diffuse  operation  is  limited  and  only  the  unstable  mode 
exists.  It  is  interesting  to  note  that  for  a  given  number  density,  the  steady  state,  stable 
surface  temperature  does  not  vary  significantly  with  the  current  density  over  a  wide 
range  (see  Fig.  3).  This  is  because  in  this  region,  the  ion  current  primarily  determines 
the  sheath  voltage  drop  which  is  then  to  a  very  good  approximation  a  constant  since 
the  number  density  is  fixed.  Heating  due  to  ion  bombardment  is  then  balanced  by  the 
external  cooling.  Therefore,  inspection  of  equation  (11)  reveals  that  the  surface 
temperature  should  not  depend  significantly  on  in  this  region.  By  contrast,  changes 
in  r.M  at  a  fixed  j„  affect  the  stable  operation  drastically  via  increased  electron  and/or 
ion  bombardment  This  is  reflected  in  Fig.  4.  The  variation  of  the  surface 
temperatures  with  electrode  work  function  is  shown  in  Fig.  5  for  Tcool  =  500  K, 
y„  =  106  A/m2,  and  n„  =  8  x  1021  m-3.  The  previously  discussed  stable  and 
unstable  regimes  are  again  found  to  exist  The  stable,  steady  state  surface  temperature 
is  almost  invariant  with  respect  to  4>  in  this  regime,  again  because  the  ion  current 
principally  determines  the  sheath  drop,  and  ion  bombardment  dominates.  Also  shown 
in  this  figure  is  the  solution  for  a  different  value  of  the  thermionic  emission  constant 
( A  =  120 x  104  Aim} IK1).  The  effect  of  a  higher  value  of  A  is  simply  to  shift  the 
turning  point  (i.e.  the  point  where  stable  and  unstable  solutions  merge  and  then 
disappear)  toward  the  higher  work  function  side.  Otherwise,  the  stable,  steady  state 
surface  temperature  is  practically  the  same  as  for  a  lower  value  of  A,  since  ion 
bombardment  is  still  the  dominant  mechanism.  Figure  6  displays  the  surface 
temperatures  versus  the  external  coolant  temperature  for  <j>  =  2.63  V,  j„  =  106  Aim2, 
n„  =  5  x  1021  m~3,  and  h  =20  kW/m2/K.  Both  stable  and  unstable  regimes  are 
seen  to  exist,  and  exhibit  the  features  that  have  just  been  discussed. 

This  section  has  focused  on  results  for  a  thermionically  emitting  cathode  operating 
under  typical  MPD  conditions.  Stable  and  unstable  regimes  as  well  as  limits  on  steady 
state  operation  have  been  found.  These  limits  are  found  to  be  strongly  influenced  by 
both  discharge  and  electrode  parameters.  The  next  section  will  focus  on  a  detailed 
discussion  of  the  operating  limits  that  have  been  found. 


4.  Limits  on  the  Steady  Diffuse  Mode 

In  this  section,  two  different  mechanisms  that  lead  to  local  melting  of  the  cathode 
surface  are  discussed.  The  first  mechanism  is  an  excessive  heating  due  to  electron 
bombardment,  that  occurs  when  the  sheath  voltage  drop  falls  below  a  critical  value.  A 
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positive  feedback  between  electron  bombrrdment  and  surface  electron  emission  results 
in  this  thermal  runaway.  The  second  mechanism  is  heating  due  to  ion  bombardment, 
which  dominates  as  the  sheath  voltage  drop  increases  sharply  for  a  relatively  small 
increase  in  the  plasma  current  density.  A  steady  state  can  no  longer  be  maintained, 
unless  plasma  discharge  parameters  change  dramatically. 


We  consider  the  thermal  runaway  mechanism  first.  The  value  of  the  critical  voltage 
drop  at  which  the  thermal  runaway  due  to  electron  bombardment  occurs,  can  be 
obtained  quite  readily  for  the  case  of  pure  thermionic  emission.  For  a  given  plasma 
current  density  and  plasma  number  density  at  the  sheath  edge,  equation  (10)  may  be 
substituted  into  equation  (11)  and  the  resulting  expression  differentiated  with  respect  to 
T  to  yield: 


dQ 

dT 


g<t>  o  ,  2kT~iE  .  etfr 
kT  eT  kT 


(13) 


For  stability  (in  the  Lyapunov  sense),  we  require  that  dQ/dT  be  less  than  zero.  A 
sufficient  condition  for  stability  from  (13)  is: 


—  >  2  (14) 

Je 

which  by  using  (7)  and  (8),  gives  a  condition  on  the  sheath  voltage  drop  Vc : 


(15) 


For  argon  and  an  electron  temperature  of  12000  K,  this  yields  an  approximate  value  of 
5.5  V  for  the  critical  sheath  voltage  drop.  It  is  important  to  note  that  this  is  only  a 
sufficient  condition  for  stability.  From  equation  (13),  it  can  be  seen  that  the  exact 
value  of  the  critical  voltage  drop  depends  on  the  electrode  properties  (work  function 
and  emission  current),  surface  temperature,  and  external  cooling  (h).  This  more 
general  criterion  on  the  total  sheath  voltage  drop  necessary  to  ensure  stability  of  a 
given  steady  state  is: 


ehT 

kjETM+etykT) 


Above  the  value  given  by  the  critical  voltage  drop,  a  steady  state  defined  by  (10)  and 
(11)  (with  Q  =  0)  can  be  maintained  indefinitely  even  for  small  fluctuations  in  n„  and 
j„.  Below  this  value,  the  positive  feedback  between  electron  bombardment  and 
surface  electron  emission  leads  to  thermal  runaway  and  a  stable  steady  state  cannot  be 
maintained.  The  surface  subsequently  melts  locally.  Alternatively,  the  condition 
dQ/dT  <  0  may  be  interpreted  as  providing  a  constraint  on  the  amount  of  cooling 
(i.e  the  value  of  h)  necessary  for  steady  and  stable  operation.  Equation  (14)  along 
with  (i9)  may  therefore  be  rewritten  in  the  following  form,  yielding  an  important 
upper  limit  on  the  local  surface  temperature  for  diffuse  operation: 


jE(T)  =  AT2exp(-e$/kT)  <  j„  -  je 


(17) 


Under  MPD  conditions,  the  critical  temperatures  predicted  by  equation  (17)  are  well 
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below  the  material  melting  temperature. 

A  second  mechanism  exists  that  can  cause  local  melting  of  the  cathode  surface.  For 
high  j «,  (given  an  nj)  or  low  n„  (given  a  /„),  it  is  possible  that  j„  — »  y,  +  jE.  When 
this  occurs,  the  total  sheath  voltage  drop  increases  sharply.  This  can  be  seen  from 
equation  (10)  where  V c  — »  ©o  as  j  ^  — >  j-t  +  jE.  This  has  two  important  consequences. 
First,  emitted  electrons  cause  increased  ionization  outside  the  sheath  since  they  gain 
large  amounts  of  energy  after  traversing  the  sheath.  Second,  the  ions  gain  substantial 
amounts  of  energy  while  falling  through  the  sheath.  These  cause  a  rapid  increase  in 
ion  bombardment  with  subsequent  heating  of  the  surface.  This  temperature  rise 
necessarily  results  in  increased  surface  electron  emission.  Consequently,  there  arises  a 
conflict  between  the  emission  current  necessary  for  cooling  the  cathode  surface  (i.e.  to 
balance  heating  by  ion  bombardment)  and  the  emission  current  necessary  to  preserve 
overall  current  conservation  at  steady  state.  The  increasing  sheath  drop  also  drastically 
reduces  the  electron  current  from  the  plasma  necessary  to  counteract  an  increasing 
emission  current,  and  hence  maintain  a  steady  state.  This  leads  to  unsteady  behaviour 
in  the  sheath,  thereby  forcing  the  surface  and  the  plasma  to  also  become  unsteady. 
The  local  surface  temperature  can  rise  due  to  ion  bombardment  in  the  meanwhile,  until 
the  surface  has  locally  melted.  Although  the  present  theory  becomes  inapplicable  at 
this  point,  it  is  capable  of  predicting  this  limit.  An  important  phenomenon  associated 
with  this  second  mechanism  but  not  with  the  first,  is  the  triggering  of  plasma 
oscillations[30].  When  the  sheath  voltage  drop  rises  and  approaches  energies  equal  to 
or  greater  than  that  of  the  first  excited  states  of  the  propellant  atoms,  ionization  levels 
can  be  enhanced  in  the  pre-sheath  region.  Also,  while  the  thermionically  emitted 
electrons  relax  their  momenta  relatively  quickly  due  to  elastic  collisions  with  neutrals, 
their  energy  relaxation  times  are  much  longer  due  to  infrequent  collisions  with  plasma 
electrons.  When  the  number  density  of  thermionically  emitted  ("beam")  electrons 
approaches  the  local  number  density  of  plasma  electrons,  longitudinal  oscillations  can 
result[30].  This  condition  may  be  quantitatively  expressed  within  the  context  of  the 
present  work  as: 

ne  <c  2.563x10 13  y/5  Vc8/5  (18) 

where  ne  is  the  number  density  of  plasma  electrons  immediately  adjacent  to  the 
sheath.  This  second  mechanism  which  is  related  to  a  high  back-EMF,  may  explain  the 
Onset  phenomenon  observed  in  the  MPD  thrusters[32,33]. 

The  two  aforementioned  mechanisms  responsible  for  limiting  the  diffuse  mode,  also 
cause  local  melting  of  the  cathode  surface.  Local  surface  melting  may  be  a  precursor 
to  spot  formation.  However,  it  must  be  pointed  out  that  an  intermediate  stable  state 
may  exist  where  molten  regions  that  are  submicron  in  size  (microspots[l])  are 
scattered  over  the  cathode  surface,  while  most  of  the  discharge  appears  diffuse.  It  is 
therefore  possible  that  the  two  mechanisms  found  and  discussed  here,  correspond  to 
transitions  from  diffuse  to  spot  mode  or  another  intermediate  mode  with  microspots. 
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5.  Effects  of  Field-Enhanced  and  T-F  Emission 

We  have  thus  far  considered  a  purely  thermionically  emitting  cathode.  In  general 
however,  the  emission  current  density  depends  both  on  the  surface  temperature  and  on 
the  surface  electric  field[31].  In  this  section,  the  influence  of  this  field-enhanced  (also 
known  as  Schottky  emission)  and  temperature-field  (T-F)  emission  is  evaluated  and 
discussed. 

The  governing  equations  for  the  sheath  and  the  cathode  surface  are  described  by  the 
surface  energy  balance  (equation  (1))  and  overall  current  conservation  (equation  (9)). 
These  equations  are  still  applicable  in  the  presence  of  an  electric  field  at  the  surface  of 
the  cathode.  However,  a  simple  analytic  solution  cannot  be  obtained  because  the 
sheath  voltage  drop  cannot  be  explicitly  determined  in  terms  of  the  surface 
temperature.  This  is  because  the  emission  current  density  in  (9)  depends  on  T  and  Ec , 
and  Ec  in  turn  depends  on  Vc.  The  required  additional  equation  for  Ec  may  be 
obtained  from  a  solution  of  Poisson’s  equation  in  the  sheath.  An  outline  of  this 
procedure  is  given  by  Prewett  and  Allen[25]  for  the  case  of  a  constant  surface 
emission  current  density.  Extension  of  their  solution  for  the  potential  distribution  in 
the  sheath  for  the  case  where  the  emission  current  density  depends  on  both  T  and  Ec , 
is  straightforward  and  will  not  be  discussed  here.  The  additional  equation  for  Ec  for 
monoenergetic  ions  is  the  following  implicit  equation: 

\  (Ec2-E„2)  =  nJT. 


2eV  ' 

1/2 

+  «•: 

-  2  +  exp[-eVclkToa] 

-  jEOmeVcle  )v2  (19) 

where  Ew  is  the  electric  field  at  the  sheath  edge  where  V  =  0.  In  general,  Ew  is  much 
smaller  than  Ec  and  may  be  ignored.  jE  is  the  emission  current  density  which  may  be 
calculated  from  quantum-mechanical  considerations.  A  general  integral  expression  for 
jE  as  a  function  of  T  and  Ec  is  given  by  Murphy  and  Good[31].  For  low  values  of 
Ec,  jE  may  be  simply  expressed  by  the  w-n  known  Schottky  formula[31].  For  the 
higher  values  of  Ec,  the  emission  is  influenced  by  both  Ec  and  7,  and  is  known  as  T- 
F  emiss'  n.  Thus,  equations  (1)  with  (2=0,  (9),  and  (17)  represent  three  coupled  and 
non-linear  equations  for  the  three  unknowns  T,  Vc,  and  Ec .  This  system  has  been 
solved  numerically.  It  is  also  possible  to  numerically  integrate  Poisson’s  equation  to 
determine  the  variation  of  the  potential  and  the  electric  field  in  the  sheath.  Figure  7 
displays  the  computed  variation  of  the  electric  field  versus  vertical  distance  from  the 
surface  for  a  particular  value  of  j„  and  n„.  However,  it  must  be  pointed  out  that  the 
resultir  ;  rofiles  represent  the  solution  within  the  sheath  only.  A  uniformly  valid 
solution  tor  the  sheath,  transition  region,  and  the  plasma  beyond  has  to  be  obtained 
through  matched  asymptotic  expansions[26,27,28].  Consequently,  the  solution  of 
Poisson’s  equation  given  here  represents  only  the  "inner"  solution.  For  the  range  of 

•7 

parameters  considered  in  this  paper,  computed  values  of  Ec  are  of  the  order  of  10 
V/m  or  less.  Such  values  of  the  electric  field  produce  a  negligible  increase  in  current 
emission  at  surface  temperatures  of  interest  to  the  MPD  thruster.  Consequently,  the 
stable  diffuse  mode  is  well  represented  bv  pure  thermionic  emission 
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6.  Discussion  and  Conclusions 

The  thermal  response  of  the  cathode  in  an  MPD  thruster  has  been  studied  under  steady 
state  diffuse  conditions.  The  electrode-adjacent  sheath  has  been  included  for 
monoenergetic  ions.  The  solution  of  Poisson’s  equation  in  the  sheath  together  with 
overall  current  conservation  and  an  energy  balance  at  the  electrode  surface,  serve  to 
determine  the  electric  field  at  the  surface,  the  surface  temperature,  and  the  total  sheath 
voltage  drop  simultaneously  and  self-consistently.  Two  operating  regimes  are  found. 
One  is  stable,  while  the  other  is  an  unstable  thermal  runaway  caused  by  excessive 
electron  bombardment.  The  thermal  runaway  leads  to  eventual  local  melting  of  the 
cathode  surface,  and  may  be  a  precursor  to  spot  formation.  The  stable  steady  state  on 
the  other  hand,  is  strongly  influenced  by  discharge  and  electrode  parameters.  The 
values  of  n„,  j„,  and  T„  used  here  have  been  obtained  from  an  approximate  non¬ 
equilibrium  model  of  the  electrode-adjacent  flow[22].  The  rate  of  cooling  is  found 
also  to  be  extremely  important.  Varying  values  of  h  result  in  two  operating  regimes, 
similar  to  the  results  displayed  in  Fig.  6.  Furthermore,  the  dependence  of  the  surface 
temperature  on  h  is  quite  significant.  For  n„  =  5  x  1021  m  "3,  7oo=106  Aim 2, 
TCOoi  =  500  K,  and  <}>  =  2.63  V,  a  1%  decrease  in  h  results  in  about  a  10%  increase  in 
the  local  surface  temperature  (see  Fig.  8).  Of  course  for  h  larger  (=  105  W lm2lK) 
than  the  value  used  here,  the  variation  is  less  dramatic  but  still  non-negligible. 

Two  mechanisms  are  found  to  be  responsible  for  destabilizing  the  stable,  steady, 
diffuse  mode.  The  first  occurs  at  values  of  the  sheath  voltage  drop  below  a  critical 
value.  Small  sheath  voltage  drops  lead  to  increased  electron  bombardment  and 
subsequent  thermal  runaway  if  the  cooling  is  insufficient.  Such  a  condition  can  occur 
in  the  MPD  thruster  operating  at  near-onset  total  currents.  Near  onset,  the  current 
density  is  sharply  peaked  and  large  near  the  inlet  and  the  exit  regions,  while  going 
through  a  minimum  in  the  middle  region  (in  the  axial  direction)  of  the  thruster.  The 
current  density  in  the  middle  region  of  the  thruster  becomes  smaller  with  increasing 
total  current  to  the  thruster,  due  to  a  high  back-EMF[32,33]  which  then  leads  to 
smaller  sheath  voltage  drops.  The  second  mechanism  due  to  excessive  ion 
bombardment  and  accompanied  by  plasma  oscillations,  occurs  at  high  values  of  the 
current  density  (for  a  given  plasma  number  density)  or  for  low  values  of  the  plasma 
number  density  (for  a  given  current  density).  Here,  high  sheath  voltage  drops  result  in 
increased  ion  bombardment  leading  to  increased  surface  electron  emission.  Since 
plasma  electrons  are  repelled  at  high  sheath  potentials,  overall  current  conservation 
given  by  (9)  cannot  be  satisfied  at  steady  state.  The  sheath  is  then  no  longer  steady. 
The  second  mechanism  may  be  inhibited  if  an  increase  in  the  sheath-edge  charged 
particle  number  density  nx  occurs  via  increased  ionization  due  to  electrons  emitted 
from  the  cathode  surface.  Clearly,  this  requires  overall  sheath  voltage  drops  of  the 
order  of  the  propellant  ionization  potential.  The  occurence  of  both  mechanisms  is 
strongly  influenced  by  discharge  parameters. 

The  simple  model  discussed  in  this  paper  represents  a  first  attempt  at  relating  the 
plasma  discharge,  the  sheath,  and  the  electrode.  Although  this  model  has  revealed 
some  important  underlying  physics,  there  are  several  deficiencies.  First,  the  pre-sheath 
has  not  been  considered  in  detail.  The  pre-sheath  is  important  not  only  for 
determining  the  ion  velocity  distribution  function,  but  also  fir  asymptotically 
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determining  a  uniformly  valid  potential  distribution.  This  potential  distribution  will 
then  yield  the  correct  variation  in  the  sheath  when  considered  on  the  scale  of  the 
Debye  length,  and  will  also  yield  the  correct  variation  in  the  plasma  when  viewed  on 
the  macroscopic  length  scale  in  the  problem.  Such  pre-sheath  solutions  exist  as 
mentioned  previously  in  section  2,  but  only  for  charge-exchange  reactions  or  in  the 
absence  of  surface  electron  emission.  For  the  plasma  in  the  MPD  thruster,  additional 
ionization  by  electrons  emitted  from  the  surface  is  important  under  some  conditions. 
Second,  radiative  heating  of  the  electrode  surface  by  the  plasma  has  been  neglected. 
This  not  only  contributes  to  direct  heating,  but  also  to  increased  heating  by  particle 
bombardment  via  increased  charged  particle  number  densities  at  the  sheath-edge. 
Radiative  transfer  in  the  MPD  plasma  is  a  topic  for  detailed  study.  Although  these 
two  deficiencies  can  affect  the  quantitative  predictions  of  the  surface  temperature,  they 
will  not  affect  the  existence  of  the  thermal  runaway  mechanism  and  unstable 
operating  regimes  found  here.  The  theory  presented  in  this  paper  together  with  models 
of  the  flowing  plasma[22,32,33],  provide  the  designer  with  a  number  of  tools.  In 
addition  to  the  prediction  of  erosion  rates  under  a  variety  of  possible  operating 
conditions,  the  proper  cathode  length,  diameter,  material,  external  cooling  conditions, 
and  choice  of  propellant  can  be  systematically  evaluated. 

Experimental  verification  of  the  predictions  of  the  theory  presented  in  this  paper 
requires  simultaneous  measurement  of  the  plasma  current  density,  plasma  charged 
particle  number  density  near  the  electrode,  amount  of  external  cooling,  and  the  local 
surface  temperature.  Although  temperature  measurements  in  low-power  steady-state 
MPD  thrusters  have  been  made[34],  lack  of  knowledge  with  regard  to  the  local  current 
and  number  densities  makes  comparison  between  this  theory  and  existing 
measurements  unreliable.  The  measured  steady  state  cathode  surface  temperatures  in  a 
subscale  device  for  thoriated  tungsten,  range  from  2087  K  to  2281  K  for  a  range  of 
total  currents  from  500  A  to  1200  A  and  mass  flows  from  46.4  mg/s  to  150.8 
mg/s[34].  This  is  comparable  to  the  parameter  range  considered  in  this  paper.  Despite 
the  difficulties  encountered  in  non-intrusively  measuring  plasma  number  densities  and 
current  densities,  efforts  to  measure  these  variables  along  with  the  local  surface 
temperature  would  be  desirable  and  invaluable.  Such  local  surface  temperature 
measurements  if  based  on  optical  pyrometry,  must  also  account  for  reflection  from  the 
radiating  MPD  plasma.  As  a  final  point,  it  must  be  mentioned  that  scanning  electon 
micrographs  (SEMs)  of  used  cathodes  always  display  evidences  of  craters  from  spots 
and  microspots  in  some  regions.  These  spots  can  easily  be  formed  during  the  start-up 
transient  or  the  shut-down  transient  of  the  arc.  Therefore,  future  experiments  aimed  at 
understanding  the  diffuse  mode  limits  must  exercise  great  caution  to  make  sure  that 
spots  are  not  formed  during  arc  initiation  or  shut-down. 
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Figure  Captions 

Fig.  1:  The  cathode  is  shown  here  along  with  the  various  heating  and  cooling 
mechanisms. 

Fig.  2:  A  typical  variation  of  the  net  heat  per  unit  area  per  unit  time  into  the  cathode 
surface  is  shown  here  versus  surface  temperature.  This  is  the  variation  predicted  by 
equations  (10;  and  (11)  in  section  2. 

Fig.  3:  The  steady  state  surface  temperature  as  predicted  by  equations  (10)  and  (11)  is 
shown  here  versus  the  net  plasma  current  density  for  various  values  of  n„.  The 

other  parameters  are  h  =  20  kWlm2/K,  <f>  =  2.63  V,  Tcool  =  500  K, 

A  =  3  x  10 4  A  lm2lK 2  and  T„  =  12000  K. 

Fig.  4:  The  steady  state  surface  temperature  as  predicted  by  equations  (10)  and  (11)  is 
shown  here  versus  the  net  plasma  current  density  n M  for  various  values  of  jM.  The 

other  parameters  are  h  =20  kW Im  2IK ,  <$>  =  2.63  V,  Tcool  =  500  K, 

A  =  3x  104  A  lm2lK2,  and  T„  =  12000  K. 

Fig.  5:  The  steady  state  surface  temperature  as  predicted  by  equations  (10)  and  (11)  is 
shown  here  versus  the  electrode  work  function  <j>  for  jeo=\06A/m2, 
n„  =  S  x  1021  m"3,  Tcool=  500  K,  h  =20  kW/m2IK,  and 

rM  =  12000  K. 

Fig.  6:  The  steady  state  surface  temperature  as  predicted  by  equations  (10)  and  (11)  is 
shown  versus  the  external  coolant  temperature  Tcool  for  <{>  =  2.63  V,  j„  =  106  Atm2, 
h  =20  kW/m2/K,  A  =  3  x  104  AIm2/K2,  and  =  12000  K. 

Fig.  7:  The  variation  of  the  electric  field  (-E)  in  the  steady,  collisionless  electrode- 
adjacent  sheath  is  displayed  here  versus  vertical  distance  from  the  cathode  surface. 
This  plot  is  a  result  of  integrating  the  Poisson  equation  with  the  limits  of  integration 
provided  by  a  simultaneous  solution  of  equations  (1)  (with  Q  =  0),  (9),  and  (17).  The 
case  displayed  here  is  for  j „=  \06  A  Im2,  n„  =  5  x  1021  m~3,  =  12000  K, 

TCooi  ~  500  K,  <{>  =  2.63  V,  h  =  20  kW/m2/K2.  The  calculated  surface  temperature  is' 
1824  K,  and  the  value  of  the  electric  field  at  the  cathode  surface  is  1.76  x  107  V/m. 
Note  that  the  electric  field  is  negative  because  of  the  convention  used  here.  The  actual 
field  is  positive,  i.e.  pointing  toward  the  surface. 

Fig.  8:  The  variation  of  the  cathode  surface  temperature  as  predicted  by  equations 
(10)  and  (11)  is  shown  here  versus  the  heat  transfer  coefficient  h,  for  a  net  current 
density  of  j„=\06A/m2,  and  for  two  different  number  densities.  The  other 
parameters  are  =  12000  K,  4>  =  2.63  V,  and  Tcool  =  500  K. 
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Nomenclature 

A  =  thermionic  emission  coefficient 
d  =  anode  thickness 

e  =  electronic  charge 

J  =  total  current 

j  =  current  density  passing  through  the  anode 

Je  =  electron  current  density  from  the  olasma  (Its  sense  is 
positive  away  from  the  anode.) 
jE  =  current  density  of  thermionically  emitted  electrons 
(Its  sense  is  positive  into  the  anode.) 
j*  -  average  net  current  density  through  the  plasma,  taken 
positive  away  from  the  anode 
k  =  Boltzmann's  constant 

L  =  anode  length 

me  =  electron  mass 

nt.  =  electron  number  density  at  the  pla.ma-sheath 
boundary 

0(  )  =  order  of  magnitude  of  the  quantity  in  parenthesis 

q0  =  heat  flux  on  the  anode  inner  surface  (Its  sense  is 
positive  into  the  anode  surface.) 

T  =  temperature 

T„  =  anode  inner  surface  temperature 

T, =  critical  value  of  Tn  where  dja,/dT0  =  0 
Td  =  anode  outer  surface  temperature 
f,  =  plasma  el jctron  temperature 

VA  =  anode  sheath  voltage  drop,  defined  as  the  potential  at 
the  plasma-sheath  edge  minus  the  anode  potential 
W  -  anode  width 

<r  =  outward  coordinate  taken  to  be  0  at  the  anode  inner 

surface,  and  d  at  the  anode  outer  surface 
(  =  '""'ssivity  of  outer  anode  surface 

X  =  thermal  conductivity  of  anode  material 
=  Stefan-Boltzmann  constant 
a  =  electrical  conductivity  of  anode  material 

<t>4  =  anode  material  work  function 


I.  introduction 


AN  important  consideration  in  ne  use  of  magnetoplasma¬ 
dynamic  (MPD)  thrusters  operating  at  steady  state  for 
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space  missions  is  electrode  erosion.  Erosion  in  MPD  thrusters 
has  been  studied  experimentally  and  has  focused  mainly  on  the 
cathode.1-3  There  are  also  ongoing  experiments  in  cathode 
erosion.4  However,  in  this  article,  we  present  a  simple  analysis 
of  the  energy  balance  on  the  anode  in  an  effort  to  predict 
anode  surface  temperatures.  This  can  subsequently  be  used  to 
predict  erosion  rates  by  evaporation.  As  will  be  shown,  even 
such  a  simple  model  reveals  subtle  physical  phenomena. 

Electrode  erosion  is  connected  with  the  modes  of  current 
conduction  through  the  electrode  surface.  Two  modes  are 
known  to  exist:  spot  and  diffuse.  The  diffuse  mode  at  subon¬ 
set  conditions  is  the  foeus  of  this  article.  Vainberg  et  al.5  have 
considered  anode  behavior  at  onset  conditions  and  beyond. 
Their  explanation  of  anode  melting  rests  on  the  anode  sheath 
reversal  mechanism  which  has  also  been  observed  by  Hugel.5 
In  constrast  to  these  earlier  works,  it  is  shown  in  this  article 
that  a  thermal  runaway  may  occur  well  before  the  sheath 
reversal.  This  work  supports  earlier  conjectures  that  local 
surface  melting  of  the  anode  precedes  spot  formation.7  Two 
operating  modes  are  predicted  by  the  theory  presented  here. 
One  of  these  yields  a  stable  steady-state  temperature  for  the 
anode,  whereas  the  other  results  in  a  thtrmal  runaway  in 
which  the  anode  regeneratively  heats  itself  until  it  melts.  The 
total  current,  anode  geometry,  and  material  work  function  are 
shown  to  strongly  influence  the  steady-state  anode  tempera¬ 
ture  for  the  stable  operating  regime. 

In  the  following  section,  the  governing  equations  describing 
anode  heat  transfer  will  be  derived.  The  solutions  to  these 
equations  under  some  conditions  of  interest  are  given  and 
discussed  in  Sec.  Ill,  followed  by  the  summary  and  conclu¬ 
sions  in  Sec.  IV. 


II.  Anode  Energy  Balance 

Consider  the  hollow  cylindrical  anode  geometry  of  the 
MPD  thruster  modeled  as  a  long  thin  slab  of  length  L,  width 
W,  and  thickness  d  (shown  >n  Fig.  1).  At  steady  state,  the 
energy  balance  gives 


(1) 


Equation  (1)  is  subjected  to  the  following  boundary  condi¬ 
tions: 


.dr. 

-x— Uo-flo 

(2) 

—  X  —  =  cosbLj 

(3) 

For  the  case  of  constant  properties,  the  system  of  Eqs.  (1-3) 
may  be  readily  integrated  to  give 


r4  ,2d 

(aSB‘d  ~ - h  do 

a 


(4) 


and 


j2d2 

X7j/  -  — — — <7(g/  +  X7o 


(5) 
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Fig.  1  The  hollow  cylindrical  anode  is  shown  in  this  idealization  as  a 
long  thin  slab  of  length  L ,  width  Wy  and  thickness  d. 


Now  ,  if  q0  is  known  andy  is  given,  Eqs.  (4)  and  (5)  represent 
two  equations  for  the  two  unknowns  T0  and  Td. 

The  heat  flux  on  the  anode  inner  surface  q0  must  be  found 
by  considering  particle  bombardment  from  the  plasma 
through  the  sheath.  Consider,  therefore,  a  power  balance  per 
unit  area  on  the  anode  inner  surface.  The  focus  here  is  on 
conditions  prior  to  sheath  reversal;  therefore,  the  anode  is  at 
a  lower  potential  with  respect  to  the  plasma  potential  at  the 
sheath  edge.  The  power  balance,  then,  mainly  consists  of  a 
balance  between  electron  bombardment  and  thermionic  emis¬ 
sion  cooling.  It  can  be  shown  that  under  these  operating 
conditions  evaporative  cooling  is  a  relatively  small  effect.  Ion 
bombardment  may  be  easily  included  in  our  analysis,  but  we 
neglect  it  for  simplicity.  Although  this  is  expected  to  influence 
final  results  quantitatively,  the  conclusions  regarding  the 
mechanisms  responsible  for  a  thermal  runaway  will  not  be 
altered.  We  then  have 

i.  =  j,  ( <t>A  +  2k  T,/e )  -  A  T0: 


200  j- 


f-  -ioo  r 
< 


ANODE  INSIDE  TEMPERATURE  (K.) 

Fig.  2  A  typical  variation  of  the  heat  defect  (defined  as  the  net  heat 
flux  into  the  anode  minus  the  net  heat  flux  out)  is  shown  here  vs  the 
anode  inside  surface  temperature  To. 

where  q0  is  given  by  Eq.  (9).  Given  d,  L,  IV,  T„  and  J,  the 
axial  current  density  throught  the  anode  ranges  up  to  j  =  J/ 
Wd,  and  the  average  current  density  through  the  plasma  at  the 
anode  surface  is  j^-J/WL.  Although  this  simplification  is 
unnecessary  (since  j  and  jx  can  be  treated  as  separate  point- 
wise  parameters),  it  allows  the  thermal  steady  state  to  be 
related  directly  to  the  total  current  (a  more  tangible  quantity) 
rather  than  the  current  densities.  With  these  quantities,  Eq. 
(10)  may  be  solved  for  T0.  Once  T0  is  determined,  Td  may  be 
readily  obtained  from  Eqs.  (4)  or  (5).  Equations  (9)  and  (10) 
have  been  solved  in  this  manner  for  representative  conditions 
in  the  MPD  thruster.  These  results  are  discussed  in  the  follow¬ 
ing  section. 


x  (c>.,  +  2kT„/e)  exp  !  ~ei>A/kT0)  (6) 

The  electron  current  density  from  the  plasma  je  can  be  deter¬ 
mined  from  overall  current  conservation 

j,  =y«  +  A  To’  exp[ -e<t>A/kT0\  (7) 

Inclusion  of  ion  bombardment  would  involve  additional  terms 
on  the  right  hand  side  of  Eqs.  (6)  and  (7).  The  anode  sheath 
drop  may  be  determined  from  Eq.  (7)  to  be 

kT 

VA~ - ~  &»[(/«  +jEVjr)  (8) 

where  jr=ent(kTt/2itme)  ■  and  jE  =  A  T,j  expj -e<t>A/kT0}. 
The  anode  is  at  a  potential  of  -  VA  with  respect  to  the  plasma- 
sheath  edge  (taken  as  the  V  -0  datum)  in  the  regime  under 
consideration.  Combining  Eqs.  (6)  and  (7)  gives 


lAkT„{Te  -  T0)  , 

+ - - - exPi 


-e<t>A/kT0\ 


(9) 

Finally,  combining  Eqs.  (4)  and  (5),  we  get  an  implicit  equa¬ 
tion  for  Tfl 


F{Tn)J-^q0-taSB(^T0-^-J^-)  =  0 


GO) 


III.  Analysis  and  Results 

In  the  previous  section,  it  was  shown  that  an  energy  balance 
on  the  anode  at  steady  state  yields  an  implicit  equation  for  the 
anode  inside  the  surface  temperature  T0.  This  equation  [Ea. 
(10))  is  nonlinear  in  T0  and  must  be  solved  numerically.  This 
section  will  focus  on  the  solution  of  this  implicit  equation, 
resulting  in  the  discovery  of  a  thermal  runaway  mode  and 
culminating  in  a  discussion  of  various  operating  limits. 

For  illustrative  purposes,  let  us  consider  a  tungsten  anode  of 
width  30  cm  and  length  10  cm.  Let  the  electron  temperature  T , 
be  fixed  at  15,000  K.  The  typical  variation  of  the  implicit 
function  F(T0)  vs  T0,  given  by  Eq.  (10),  is  displayed  in  Fig.  2. 
The  function  F(T0)  represents  the  net  heat  into  the  anode  per 
unit  area  per  unit  time.  It  is  immediately  evident  from  Fig.  2 
that  two  steady-state  solutions  exist.  Without  performing  a 
stability  analysis,  one  can  conclude  on  the  basis  of  physical 
reasoning  that  the  lower  temperature  is  a  stable  root  and  the 
higher  one  is  an  unstable  root.  Consider  the  smaller  tempera¬ 
ture  first.  Any  perturbation  to  the  left  of  this  root  (i.e.,  on  the 
lower  temperature  side)  results  in  a  net  heat  flux  into  the 
anode.  To  counter  this,  the  surface  temperature  must  increase 
in  order  to  maintain  a  steady  state.  Similarly,  any  small  per¬ 
turbation  to  the  right  (i.e.,  on  the  higher  temperature  side) 
results  in  a  net  heat  flux  out  of  the  anode.  To  counter  this,  the 
surface  temperature  must  decrease.  Thus,  we  see  that  the 
tendency  of  the  thermal  response  is  to  drive  the  temperature 
toward  this  root.  Using  this  physical  argument  on  the  second 
root  (i.e.,  higher  temperature),  we  may  conclude  that  it  is 
unstable.  If  the  surface  temperature  increases  just  above  the 
value  given  by  the  second  root,  the  steady-state  energy  balance 
indicates  a  further  increase  in  temperature  in  order  to  main¬ 
tain  a  steady  state.  This  is  because  the  only  possible  way  of 
losing  energy  is  through  radiation  and  by  thermionic  emission. 
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both  of  which  increase  with  increasing  temperature.  This  re¬ 
sults  in  a  thermal  runaway  with  eventual  melting  of  the  anode. 

A  physical  interpretation  of  the  two  aforementioned  solu¬ 
tions  will  now  be  discussed.  Examination  of  Eq.  (8)  indicates 
that  the  anode  sheath  drop  decreases  as  the  net  plasma  current 
density  increases.  This  results  in  increased  electron  bombard¬ 
ment,  as  can  be  seen  from  Eqs.  (6)  and  (7),  causing  heating  of 
the  anode  surface.  Two  possible  scenarios  may  ensue.  The 
resulting  increase  in  the  surface  temperature  leads  to  an  in¬ 
crease  in  the  emission  current  density,  which  further  lowers 
the  anode  sheath  voltage  drop.  This  positive  feedback  may 
repeat  itself  until  thermal  runaway  leads  to  local  melting  of  the 
surface.  The  second  scenario  is  stable,  steady  operation  arising 
from  the  fact  that  sufficient  cooling  of  the  anode  prevents  the 
thermal  runaway  by  excessive  electron  bombardment.  The  two 
solutions  found  to  Eq.  (10)  represent  these  two  situations. 

The  anode  surface  temperature  T0  for  stable  operation  is 
found  to  be  strongly  dependent  on  the  total  current,  anode 
geometry,  and  the  anode  work  function.  Figure  3  shows  the 
anode  thickness  for  various  total  currents,  and  Fig.  4  shows 
the  anode  inside  surface  temperature  vs  total  current  for  vari¬ 
ous  anode  thicknesses.  From  these  results,  it  is  evident  that  for 
a  given  total  current,  there  exists  an  optimum  thickness  for 
achieving  a  minimum  anode  temperature.  An  interesting  fea¬ 
ture  of  the  pair  of  solutions  found  to  the  steady-state  heat 
balance  equation  is  that  the  stable  and  unstable  roots  start 
moving  toward  each  other  as  the  current  is  increased.  This  con¬ 
tinues  until  a  critical  value  of  the  current  is  reached,  beyond 
which  no  steady  solution  can  be  found  below  the  melting  point 
of  the  material.  Also,  for  a  fixed  total  current,  these  roots 
approach  each  other  as  the  thickness  is  increased.  It  is  possible 
that  the  analytical  behavior  of  this  unstable  thermal  runaway 
point  will  yield  insight  into  electrode  material  breakdown. 

Some  interesting  limits  may  be  readily  obtained  by  analyz¬ 
ing  the  equations  presented  in  Sec.  II.  For  instance,  from  Eq. 
(5)  it  is  clear  that 

j'-d'-/2a<kT0  or  j  <jc  =  (2koT0)  ’/d  (11) 


prevent  thermal  runaway.  This  thermal  runaway  can  occur 
under  the  steady,  diffuse  mode  operation  of  the  anode. 

An  important  stability  criterion  may  be  derived  by  requiring 
that  8F/5T0<  0  for  stable  diffuse  operation.  Using  Eq.  7,  we 
find 


t-T 

Va*— —U 
e 


y»  +  0 


2  O-i  Te 


(12) 


For  VA  below  the  value  given  by  the  right-hand  side  of  Eq. 
(12),  excessive  electron  bombardment  will  cause  a  thermal 
runaway.  Thus,  local  surface  melting  can  occur  prior  to  sheath 
reversal. 

A  limit  may  also  be  found  for  the  plasma  current  density./',. 
Considering  thin  anodes  (i.e. ,  T=T0  throughout  the  anode), 
neglecting  ohmic  heating  and  any  external  cooling  except  for 
radiation,  we  may  write  the  energy  balance  at  steady  state  as 

tonTS  -  To)  exp(-eoA/kT0) 


(13) 


Using  the  fact  that  Te>T0,  differentiating  with  respect  to  T0, 
and  setting  dj^/dTQ  =  0,  we  obtain 


[1  —(4kToc/e<t>A)] 
U  [4>A+(2kTe/e)) 


eoshToc* 


(14) 


where  (/«,) is  the  extremum  plasma  current  density,  and 
is  the  critical  value  of  T0  where  dja/dT0  =  0.  Additional  heat 
transfer  to  an  external  coolant  results  only  in  a  slight  modifi¬ 
cation  of  Eq.  (14): 


[1  —  (4kToc/e<t>A)\  4  (1  -(kT^/eO.,)} 
=  ld>A+(2kTt/e)}  e<7sB  0c  [«,  +  (2kT,,/e)\ 


hToc 


which  is  very  similar  (within  a  proportionality  constant)  to  the 
thermal  runaway  condition  obtained  by  Hantzsche*  for 
cathode  spots.  However,  it  must  be  pointed  out  that  Eq.  (11) 
is  obtained  here  for  an  anode  of  constant  electrical  conductiv¬ 
ity,  operating  in  the  diffuse  mode.  This  sharply  differs  from 
Hantzsche’s  consideration  of  a  cathode  spot  with  an  electrical 
conductivity  that  varies  according  to  the  Wiedemann-Franz 
law  (i.e.,  <joc  1/7").  Furthermore,  in  constrast  to  Hantzsche’s 
work,  surface  cooling  due  to  thermionic  emission  does  not 


(15) 

where  h  is  the  heat  transfer  coefficient  between  the  anode  (at 
T0)  and  an  external  coolant  (at  T,  <*  To).  For  Te  =20,000  K, 
Eq.  (14)  approximately  yields  47.6  A/cm:  for  a  tungsten  an¬ 
ode  (<t>A  =4.52  V)  and  5.4  A/cm:  for  a  thoriated  tungsten 
anode  (<t>A  =  2.63  K).  The  lower  work  function  yields  a  lower 
maximum  plasma  current  density  because  of  a  lower  value  of 
Toe-  Since  these  limiting  current  densities  are  far  lower  than 
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Fig.  3  The  anode  inside  surface  temperature  is  plotted  here  vs  anode 
thickness,  for  various  total  currents. 
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Fig.  4  The  anode  inside  surface  temperature  is  shown  here  vs  total 
current,  for  various  anode  thicknesses. 
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those  required  for  a  steady  MPD  discharge,  it  is  clear  that  the 
anode  must  be  cooled  externally.  This  is  done  in  practice.1-4 

IV.  Summary  and  Conclusions 

An  anode  energy  balance  has  been  performed  that  includes 
the  effects  of  the  anode  sheath.  Results  show  that  under  many 
conditions,  two  steady-state  solutions  can  be  found.  One  of 
these  corresponds  to  a  stable  operating  point,  the  other  to  a 
thermal  runaway.  The  stabu  ’•c.ot  that  gives  the  anode  inside 
surface  temperature  at  steady  state  was  found  to  be  strongly 
dependent  on  the  total  current,  anode  geometry,  and  the  mate¬ 
rial  work  function.  A  stability  condition  in  the  form  of  a 
minimum  anode  sheath  voltage  drop  [Eq.  (12)]  has  also  been 
given. 

Several  conclusions  can  be  derived  from  the  anode  thermal 
analysis: 

1)  There  exists,  for  a  given  current  density,  an  optimum 
thickness  for  achieving  a  minimum  anode  temperature. 

2)  Under  many  conditions,  there  exists  a  pair  of  steady-state 
solutions  to  the  heat  balance,  one  of  which  is  a  stable  operat¬ 
ing  point,  the  other  of  which  is  an  unstable  thermal  runaway 
point. 

3)  There  exists  a  maximum  steady  operating  current  density 
for  a  given  electrode  geometry  based  on  melting  temperature 
and  cooling  considerations. 

The  thermal  runaway  mode  discovered  for  the  anode  may 
explain  anode  spots  under  MPD  conditions.  Although  some 
limits  have  been  found,  anode  heat  transfer  effects  need  to  be 
explored  further  in  order  to  learn  more  about  what  controls 
the  limits  of  stable  operation.  Since  external  cooling  is  seen  to 
be  extremely  important,  detailed  experimental  results  quan¬ 


tifying  cooling  rates  for  various  operating  conditions  and  ge¬ 
ometries  would  be  valuable.  Radiative  heat  transfer  from  the 
plasma,  which  has  been  neglected  here,  is  another  important 
area  for  further  research. 
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